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Abstract. S-type stars are late-type giants whose atmosphere is enriched in carbon 
and s-process elements because of either extrinsic pollution by a binary companion or 
intrinsic nucleosynthesis and dredge-up on the thermally-pulsing AGB. A large grid of 
S-star model atmospheres has been computed covering the range 2700 < r e ff(K) < 
4000 with 0.5 < C/O < 0.99. ZrO and TiO band strength indices as well as VJHKL 
photometry are needed to disentangle T e jf, C/O and [s/Fe]. A "best-model finding tool" 
was developed using a set of well-chosen indices and checked against photometry as 
well as low- and high-resolution spectroscopy. It is found that applying M-star model 
atmospheres (i.e., with a solar C/O ratio) to S stars can lead to errors on T e ff up to 400 K. 
We constrain the parameter space occupied by S stars of the vast sample of Henize stars 
in terms of T e ff, [C/O] and [s/Fe]. 



1. Introduction 



The S class was oris inally defined bv lMerrilll dl922l) to designate a group of curious red 



stars which did not fi t well into either class M (TiO stars) or classes R and N (carbon 
stars). iKeenanl ( 19541) clarified the situation by accepting as S stars only those exhibit- 



ing ZrO bands. The numerous attempts to link phenomenological spe ctral classifica 



tion criteria to physical parameters (T F ff, gravity, C/O, ls/Fe], [Fe/H]) dKeenan|[l954l : 
Keenan & McNeill 1 1974 lAkdl 19791 : iKeenan & Boeshaar| [l980) only lead to imprecise 



results, because low -resolution diagnostics are strongly entangled in terms of T e ff, C/O 
and [s/Fe] variations. The o nly in-depth dis cussion of the thermal structure dates back 
to the pioneering paper of iPicciriliol (Il980l) . He already insisted on the strong influ- 
ence of the C/O ratio on the atmospheric structure and spectra of S stars, in addition 
to effects due to the s-process elements overabundance. His investigation was however 
mostly limited to qualitative statements, due to obvious technical limitations. Most 
subsequent analysis of S stars relied on models designed for M-type stars, not allowing 
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for C/0 or [s/Fe] ratio changes. In the present paper we presen t a new grid of model 
atmospheres, superseeding the one presented in lPlez et"al1 (120031) . covering most of the 
parameter space of S-type stars, and attempt to provide a calibration of photometric 
indices in terms of T e ff, C/O and [s/Fe]. 



2. Model atmospheres and spectra 
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Figure 1. Thermal and pressure structure of models with T e ff = 3000 K and 
[s/Fe]=+2 dex 



Since models for S star atmospheres are vi rtually non-existent, a grid of MARCS 
model atmospheres (see lGustafsson et all ( 2008h for details on the models computation) 
for S stars has been calculated: 2700 < T eff (K) < 4000 (step of 100K); C/O = 0.5, 
0.750, 0.899, 0.925, 0.951, 0.971, 0.991; [s/Fe] = 0., +1., +2. dex ; [Fe/H] = -0.5 and 
0. dex ; log(g) = 0,1,2,3,4,5. 

All models were computed for M = 1 M Q and with [a/Fe] - -0.4 x [Fe/H]. 
Opacities as complete and accurate as possible were included, including polyatomic 
molecules and a specific ZrO linelist (described in Plez et al., in preparation). Models 
were computed through opacity sampling with more than 10 5 wavelength points, local 
thermodynamic equilibrium, mixing-length theory of convection and spherical symme- 
try for log(g) < 2. 

A total of 3522 converged model atmospheres was obtained. The model structure 
for T e ff =3000 K and [s/Fe] =+2 dex models is shown in Fig. 1, where the major influ- 
ence of C/O on the thermal structure is readily apparent (whereas the [s/Fe] ratio has 
less importance). The P gas - T500 relation (fixed mostly by log g) stays basically un- 
changed, whereas the T - tr oss relation (governed by the energy balance requirement) 
reaches higher temperatures at the surface for higher C/O. When C/O increases, P gas at 
a given T increases. The latter effects are due to a large decrease of the partial pressures 
of H2O and TiO, two major opacity contributors. Fig. |2]illustrates how the depth of TiO 
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Figure 2. Top: absorption spectra of TiO and ZrO at 3000 K. Bottom: Thick 
lines: spectra with increasing C/O ratios (from bottom to top) for T e ff=3000 K and 
[s/Fe]=+l dex. Thin dashed line: spectrum with C/O=0.99 but no s-process en- 
hancement ([s/Fe]=0 dex) 



and ZrO bands decreases with increasing C/O, for different models of T e ff=3000 K, and 
the influence on the spectra of the level of s-process enhancement. 



3. Confronting the models with observed color and spectral band indices 



Synthetic spectra are now compared to obs ervations. T he Henize sample of S stars 
(205 stars with R < 10.5 and 5 < -25°: lHenizelll960l) is of particular interest for 
that purpose, since (i) it co llects S stars with no bias against high galactic latitudes 



( Van Eck & Jorissen 2000bl). and (ii) a large observational material has been collected 
for this sample dVan Eck et al.ll2000r) . Fr om these data, t he (V - K)o, (J - K)o color- 



color diagram, dereddened according to iDrimmel et al.l (|2003), has been constructed 
(Fig. [3]>. Similarly, a set of TiO and ZrO band-strength indices have been computed 
from the low-resolution spectra, and displayed on Fig. [3] Their comparison with the 
model values make it possible to estimate T e ff , C/O and [s/Fe] since: (i) the (V - K)o, 
(J - K)o color-color diagram disentangles T e ff and C/O; (ii) the (TiO, ZrO) diagram 
disentangles T e ff and [s/Fe]. In both cases, there is a good segregation between M and 
S stars with, however, some degeneracy between C/O and [s/Fe], especially for low 

T e ff- 

The (V-K)q, (J-K)o color-color diagram reveals that, for a given V—K, the range 
in T e ff covered by models of different C/O ratios can be as large as 400 K. Therefore, 
the application to S stars of the usual M-star temperature scale based on the V—K index 
(as done in the past when specific S-star models were unavailable) leads to errors on 
Teff of up to 400 K. 
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Figure 3. Top panels: Comparison between color indices of observed M (squares), 
S (triangles) and C (crosses) stars, and color indices computed on synthetic spectra 
of S stars for [s/Fe] = (left) and 2 (right). The models with the lowest tempera- 
ture (2700K) and highest C/O ratio (0.991) are on the top of each "grid". Bottom 
panels: ZrO index versus TiO index for [s/Fe] = (left) and 1 (right). The grid cor- 
responds to solar-metallicity, log g=0 models ranging from T e ff=4000K, C/O=0.5 
(around coordinates 0.05, 0.05 on the leftmost figure) to T eff =2700K, C/O=0.99 
(around 0.3,0.65). Stars clumping around (TiO,ZrO) = (0,0) are G and K giants. 
All S stars to the left of the region covered by the grid are SC stars 



4. The atmospheric parameters of S stars 



We have built a "best model finding tool", thanks to an appropriate weighting of well- 
chosen photometric and narrow-band indices (dereddened Geneva and VJHKL photom- 
etry, ZrO, TiO and NaD band strengths) and chi-square minimization between observed 
and synthetic indices. The adequacy of the selected mode ls has been checked on low- 
resolution spectra, dereddened according to lCardelli et al.l(ll989l) : the agreement is very 
good in most cases (see Neyskens et al., this volume). 

Fig. |4] presents the distribution of Henize S giants in terms of temperature and 
C/O ratio. The temperature difference between Tc-poor (polluted binary) S stars and 
the cooler Tc-rich (genuine TPAGB) S stars is clearly visible. Among Tc-rich stars, 
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Figure 4. Comparison of the T e g distribution of Tc-rich (shaded histogram) and 
extrinsic (unshaded histogram). The 7 top panels separate the stars according to the 
C/O ratio 



despite the small number statistics, the expected gradual increase of the C/O ratio as 
the star cools down and ascends the TPAGB is also visible. 

This new grid of model atmospheres is an unavoidable prerequisite to reliable 
spectroscopic chemical analyses of these objects enriched in s-process nucleosynthesis 
products. It will allow us to pursue on a more quantitative basi s the co mparison between 
extrinsic and intrinsic S stars initiated by I Van Eck & Jorissenl (|2000al) . 
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